Improved phase gate reliability in systems with neutral Ising anyons 
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Recent proposals using heterostructures of superconducting and either topologically insulating or semicon- 
ducting layers have been put forth as possible platforms for topological quantum computation. These systems 
are predicted to contain Ising anyons and share the feature of having only neutral edge excitations. In this note, 
we show that these proposals can be combined with the recently proposed "sack geometry" for implementation 
of a phase gate in order to conduct robust universal quantum computation. In addition, we propose a general 
method for adjusting edge tunneling rates in such systems, which is necessary for the control of interferometric 
devices. The error rate for the phase gate in neutral Ising systems is parametric ally smaller than for a similar 
geometry in which the edge modes carry charge: it goes as T A rather than T at low temperatures. At zero 
temperature, the phase variance becomes constant at long times rather than carrying a logarithmic divergence. 



A device implementing topologically protected quantum 
computation would allow the fault tolerant manipulation of 
quantum information stored in the non-local properties of 
non-Abelian anyons. 12 Unfortunately, such a device is dif- 
ficult to create, as the non-Abelian anyons that are most 
likely to be accessible in a laboratory setting, Ising anyons, 
are not associated with topologically protected operations 
that are universal for quantum computation. However, Ising 
anyons can store quantum information, and their manipula- 
tion can produce a protected set of operations (the Clifford 
gates) that are sufficient for any classical operation^ In or- 
der to make the jump to quantum computation, a single qubit 
phase gate (or a gate equivalent through Clifford operations) 
is necessary. 6 Fortunately, by using a so-called "magic-state 
distillation" procedure we can tolerate a high degree of er- 
ror (14%) in the phase gate if the other basic operations are 
protected^ 

Recently, we proposed a method for implementing a single- 
qubit phase gate in an Ising-anyon-based quantum computa- 
tional platform. 7 This was the first proposal that could rea- 
sonably meet the error threshold set forth by Brayvi 5 in the 
context of this architecture. In principle, using this gate, a 
v = 5/2 fractional quantum Hall (FQH) state could be made 
universal for quantum computation. However, building quan- 
tum devices in v = 5/2 FQH systems presents several prob- 
lems. First, the v = 5/2 state itself is not fully understood, 
although recent evidence seems to point to its non-Abelian 
nature. 810 Second, high quality samples are necessary for the 
observation of the plateau at v = 5/2, as well as ultra-low- 
temperature environments; the proposed phase gate, for in- 
stance, has an error rate that grows as T at low temperatures. 
Finally, stray quasiparticles are likely to exist in this system 
even under ideal experimental conditions. In fact, their ex- 
istence is the basis for the experiment purported to show the 
non-Abelian nature of excitations at v = 5/24^ In the context 
of quantum computing, these stray quasiparticles would inter- 
fere with the fine control of the system necessary to complete 
a full quantum algorithm. 

New heterostructure systems based on s— wave supercon- 
ductors (SC) sandwiched with either topological insulators 
(TI) or semiconductors with strong spin-orbit interactions 
have been proposed recently; the excitations in these sys- 
tems are predicted to be Ising anyonsjii— These systems may 



have certain advantages compared to the other candidates for 
hosting non-Abelian excitations, in particular the v = 5/2 
fractional quantum Hall liquid, because they involve mate- 
rials that may not have to meet requirements as stringent as 
those necessary for the FQHE. One characteristic of all these 
systems is that they have chiral edge modes. However, un- 
like in FQH systems, the edge modes with non-Abelian any- 
onic nature in these new systems are neutral. This presents 
a challenge because the interferometry experiments purported 
to show the non-Abelian nature of the anyonic excitations in 
the v — 5/2 stated 9 rely on the electric signal due to their 
charge. However, recent theoretical developments have shown 
that, given the proper geometry, the interference of the neu- 
tral edge modes in the heterostructure systems may be electri- 
cally detected as well.— ^ This makes relevant the question 
of the practicality of such systems for quantum information 
processing, and in particular warrants a reinvestigation of the 
phase gate for neutral systems. We find that in systems that 
have only neutral edge excitations the error rate of the phase 
gate actually improves significantly. However, these systems 
present a separate challenge in that the edge of the system 
is not so easily deformed as it is in a FQH system; our ear- 
lier proposal for the phase gate relied on being able to shape 
the edge by electrostatic gating in order to control tunneling 
rates. Here, we propose an alternative method of controlling 
such tunneling. 

We shall begin by reviewing the analysis of Ref. [7] in the 
context of a system with a single neutral edge mode. Fol- 
lowing that proposal, we consider a topological qubit encoded 
in a pair of anyons, each carrying Ising topological charge 
a whose two possible fusion channels / and rp form a com- 
putational basis. In superconducting heterostructure systems, 
the anyons comprising the qubits may be localized using su- 
perconducting vortices (or non-superconducting regions pen- 
etrated by half a quantum of fundamental magnetic flux hc/e). 
These two anyons are placed in the "sack" geometry as shown 
in Fig.[T] 26 The quasiparticles moving along the edge can tun- 
nel from —a/2 to a/2 (so the sack has perimeter length a). 
Generally, there could be more than one type of excitations 
tunneling between —a/2 and a/2. In the weak-tunneling, low 
temperature regime, the quasiparticles with the most relevant 
tunneling operators will dominate the tunneling current. How- 
ever, quasiparticles that do not carry a will have no effect on 
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FIG. 1 : A top view of the proposed phase gate device in heterostruc- 
ture systems with Ising anyons. A hatched region corresponds to the 
region coated with a superconductor; the edge state is formed at the 
boundary of this coating with a Zeeman-inducing magnetic material. 
The SC-coated sack like region is pierced by two superconducting 
vortices, each binding a a anyon. A crosshatched region depicts a 
weak Josephson link, which is biased by the supercurrent / in order 
to control the quantum phase slip rate (see text for details). 



the topological qubit here and therefore will not enter in our 
analysis. 

Let us pause here in order to clarify an important point 
which appears to be somewhat obscure in the context of su- 
perconducting heterostructure systems. The chiral Majorana 
fermions that are present on the edge of these systems are not 
non-Abelian excitation, while Majorana zero modes bound to 
the vortices are. In other words, it is the presence of vor- 
ticity, or a twist in the boundary conditions for the Majorana 
fermions, which is essential here. It is this twist that is en- 
coded in the field a. Therefore, the excitations that must tun- 
nel in our device in order to effect the phase gate are vortices. 
This point has been made in Refs. Il5lll7ll . albeit it was orig- 
inally expected that these would be Abrikosov vortices. The 
serious drawback of Abrikosov vortices is that they tend to 
behave as very classical objects, which makes the prospects 
of their quantum tunneling or interference rather problematic. 
To the best of our knowledge, such quantum effect have not 
been observed in any systems. On the other hand, Josephson 
vortices can be quantum: their tunneling is known as quantum 
phase slips and their quantum interference has been observed 
in the context of Aharonov-Casher effect42i2£ Therefore the 
device should be fabricated with a Josephson weak link be- 
tween —a/2 and a/ 2; we should rely on controlled tunneling 
of Josephson vortices along this link. We will revisit the is- 
sue of manipulating their tunneling amplitude later, for now 
we shall return to our analysis and show how the interference 
between the possible trajectories from left to right enacts a 
non-trivial transformation on the qubit. 

The combined edge and qubit system is described by the 
Hamiltonian 



this is the Aharonov-Casher phase induced by taking a vortex 
around a region with odd electron parity^ 

The density matrix of the combined system is represented 
by x an d the qubit's density matrix is obtained from this by 
tracing out the edge p = TveX- We assume that the edge and 
the qubit are initially unentangled. Closely following Ref. 0], 
we solve the interaction picture Schrodinger equation 



i^ = [H tnn {t)®* z ,x{t)] 



where A(t) 
P(t) 



Et A{t)e~ iHEt , to obtain 
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Here 9 is the relative phase accumulated between the two basis 
states, while q 2 represents the loss of coherence due to vari- 
ance in the phase. The diagonal elements of the qubit density 
matrix are unaltered from their initial state, as the Hamiltonian 
commutes with t <g> a z . 

Computing the values of 9 and q 2 to second order in the 
tunneling Hamiltonian, we have 



2 J dt' (iftun(O) , 
Q2 



(4) 



4j dhj dt 2 {H tnn (t x )H tnn (t 2 )). (5) 

To compute concrete values of 9 and <, 2 , we use the field 
theoretic description of the edge of a system with a neutral 
Ising chiral edge mode. The Lagrangian for the unperturbed 
edge is 



(6) 



L E = i dxip(d t + vd x )tp, 



where ip describes the chiral fermion mode. The operator that 
tunnels a quasiparticles across the constriction is 



H tuXL = Te-^ +h.c. 



(7) 



where /3 includes the dynamical phase acquired in traveling 
around the sack as well as any Abelian braiding statistics fac- 
tors. 

Assuming the edge was initially in thermal equilibrium (i.e. 

X(0) = — e r _H g /r l ® /°(0)) a t temperature T, we find 



tftun(t) ) = 2 



XwT/v 
sinh^ 



1/8 



|r|cosv?, 



(8) 



H = H E ®t+H tmi (t)®tr z 



(1) 



where He is the Hamiltonian describing the unperturbed edge 
and H tun describes tunneling of Josephson vortices across the 
constriction. Here the <r z represents the braiding statistics of 
the associated a with the qubit, picking up a minus sign each 
time the a braids around the ip charge. In SC heterostructures 



Here A is a short range cutoff and <p — arg {V} — (3. 

In order to estimate <r 2 , we assume T is abruptly turned on 
to some constant value for the duration of time from to t, 
and all other quantities are held fixed. In this case, we have 

e = u 2 [ dtj dt 2 r)(t 2 - ti), (9) 
Jo Jo 
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where u = 2 ( H tU n ) while the tunneling is on, and 



with 



= J T(f) -l/4 + T(f) l/4 _ ^ 



sinh 2 (2?21) 

= i - : "/ 



sinh 2 ((i<5 - f)7rT) 



(10) 



(ID 



Note that rj(T,t) — s> exponentially for long times at a 
rate proportional to the temperature. Therefore, if the gate is 
performed slowly in comparison with the sack time a/v then 
the variance grows linearly: 

/OO POO 
At'uj 2 ri(t') - / At' \t'\u} 2 r]{t'), (12) 
-oo J — oo 

More precisely, we have for t > a/v 

e = uj 2 J At 1 (t- \t'\)r)(t'), 

/OO /» OO 

dt' (t - |i'|) 77(7/) + 2w 2 / dt' (t' - i) r/(i') 
-OO J t 

= \t + <; 2 + E(t), (13) 

where E(t) = 2lo 2 ff°At' (f ~ t) 77(f) tends to at long 
times. 

For t a/v, we may approximate 



r)(t,T)< v (t,0) 



O 



\vt 



so that 



(14) 



(15) 



This confirms that for long times t S> a/v the phase variance 
approaches the linear form ( TTZi i. At temperature, the ap- 
proach is a power law going as t~ 2 . At finite temperature, the 
approach will be even faster, with the nonlinear term decreas- 
ing exponentially. 

The long-time decay rate A when only the neutral mode is 
present is given by 



A„ = uj 2 \ At'r](t'). 



(16) 



At T = 0, this integral is exactly 0. This may be seen by not- 
ing that the integrand 77 has no poles in the lower half of the 
complex t plane when T = 0, and falls off as |f | 4 for large 
\t\. To calculate A for non-zero T, we first make a change of 
variables u = sinh(7rf'T), so that in the lower half of the com- 
plex u plane 77 has only one branch cut running from u = —i 
to u = — ioo and no poles. This allows us to change the inte- 
gral from one running along the real axis to one running along 



the branch cut. If T <C v/ira, the integral may be approxi- 
mated as 



w 2 f°° , 2 
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Where T n = 7;/7ra is the characteristic temperature scale for 
the neutral mode in the sack. This may be compared with 
the formula for a general Ising edge with charge and neutral 
modes 



A = A„ H — 

TFT 



T 2 1 (T 
t= tan ip + - — 
T 3 VT 



oil 



(18) 



when t ^> 1/ttT 3> a/v c , n , and where T = w/7ra 
(1/-D 2 = (g n - 1/8)/tj 2 + g c /v 2 ), with u Ci „ and 5 Cj „ the ve- 
locities and scaling dimensions for the charge and neutral 
modes, respectively. This formula for A may be derived us- 
ing the above method with the more general form of 77 derived 
in Ref 17[]. Note that when only the Ising neutral mode is 
present, g n = 1/8 and g c = 0, so T = 00. This eliminates 
the first order term in T and the explicit dependence on the 
Aharonov-Bohm phase <p. 

We can see that having only a single Ising neutral mode 
provides a significant advantage over having charge and neu- 
tral modes both present. First, there is a marked decrease in 
the decay rate for neutral-mode-only systems, with A ~ T 3 
rather than T at low temperatures. Second, the error does not 
depend explicitly on the Aharonov-Bohm phase, so there is no 
runaway error as the phase advancement slows near ip = ir/2 
as there is in the case where charge modes are present. Finally, 
the zero-temperature limit of the neutral system exhibits only 
a constant level of degradation in the qubit state at long times, 
while in general one would expect the off diagonal elements of 
the density matrix to decay as a power law. That is, in general 
H(i) + <Tg has a logarithmic divergence for long times when 
T = 0. In a neutral system, H(t) + ?q converges as t — » 00. 

The above analysis suggests that the most important tem- 
perature scale for a device intended to implement quantum op- 
erations via interferometry on a system with neutral Ising edge 
modes is T n = hv/kBair, where a is the typical perimeter size 
of a device element such as a sacki2£ or an interferometeriiSiS 
The edge mode velocity v is determined by the Fermi velocity 
in the host material (for devices based on topological insula- 
tors) or by the Rashba spin-orbit coupling constant (for semi- 
conductor devices), as well as by how closely the chemical 
potential is tuned to the Dirac point. 

Fu and Kane 16 calculated the edge mode velocity in a SC/TI 
heterostructure to be 



v F 



yfi -/i 2 /M 2 
1+/7.2/A 2 



(19) 



where M is the energy of Zeeman splitting and A is the in- 
duced pairing potential due to the superconducting proximity 
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effect. In order to derive an upper bound on the operating tem- 
peratures of such systems, we assume that \i has been tuned 
to near zero and that the typical device perimeter a is of order 
1 /jm. Then for the SC/TI system Be2Se3, which has Fermi 
velocity 4.6 x 10 5 m/s, the maximum operating temperature 
T n < 1 K. For the proposed SC/semiconductor system, we 
evaluate the edge mode group velocity given in Ref. lfl8ll . i.e 
v — a J dx (4>\ a x r z \<f>), where \4>) is the Nambu spinor solu- 
tion of the Bogoliubov-de Gennes equations for the Majorana 
mode. For InAs, Ref. 1 12] suggests the values of .5 meV for 
effective pairing potential and 1 meV for the Zeeman split- 
ting energy. Using these parameters and the measured Rashba 
coupling 1.4 x 10 4 m/s and effective mass m* ~ .05m e ri 
we have calculated the value of the integral to be .65. This 
leads to a velocity of 9 x 10 3 m/s. The maximum operat- 
ing temperature is therefore T n < 20 mK. In this sense the 
SC/TI system is the more promising proposal, though it may 
be possible to significantly increase the operating tempera- 
ture in the semiconductor case by using a system with a giant 
Rashba coupling. (A Bi/Ag(001) interface, for instance, has 
been measured to have a Rashba coupling of 1.2 x 10 6 m/s.— 
This corresponds to T n < 3 K.) 

A new challenge presented by the proposed systems with 
neutral Ising anyons is that, unlike quantum Hall edges, their 
edges cannot be manipulated electrostatically by gating. E.g., 
in the case of TI based systems, the edge is defined by the 
boundary between two distinct coatings: those by a supercon- 
ductor and a ferromagnet. It is extremely difficult to imagine 
how one would move such a boundary "on the go" in order 
to manipulate the tunneling rate of a particles (in this case, 
Josephson vortices) - a key to implementing our proposal for 
a phase gate. A use of multiferroics might be of some utility 
here, yet, to the best of our knowledge, no device with a tune- 
able FM-SC boundary has been fabricated to the date. What 
we envision instead is using a relatively wide Josephson junc- 
tion (to suppress the "ambient" vortex tunneling) and biasing 
it by current in order to manipulate the shape of the tunneling 
barrier. Indeed, the quantum dynamics of a Josephson junc- 
tion is usually modeled in terms of a motion of a particle in a 
"washboard" potential U(<f>) = — [/o(2/0// c + cos2</>) where 



Uq = HI c /2e with I c being the critical current. The current 
across the junction provides the tilt which in turn makes the 
potential barrier for a phase slip both lower and narrower. In 
the language of Josephson vortices the same phenomenon is 
described in terms of a Magnus force acting on a vortex due 
to the transverse current; the effect of such a force is to tilt the 
potential barrier for the vortex. Since the tunneling rate de- 
pends exponentially on the parameters of this barrier, current 
biasing should provide a simple and reliable way of manipu- 
lating such a rate for the purpose of effecting the phase gate. 

To conclude, we have discussed both advantages and disad- 
vantages arising from using the sack geometry for effecting a 
phase gate in the systems with neutral Ising anyons such as re- 
cently proposed heterostructure systems where either a topo- 
logical insulator or a semiconductor with a strong spin-orbit 
coupling (and an additional source of Zeeman splitting) is in- 
terfaced with a superconductor. One significant advantage is 
a parametrically lower error rate while the most apparent dis- 
advantage is a lack of developed techniques of manipulating 
the shape of the edges. The latter was the key to affecting the 
tunneling rates across the constriction in the original proposal. 
In this note we have proposed an alternative way of manipu- 
lating the tunneling rate using a weak Josephson link biased 
by current. 

We should also note that while this paper was in prepara- 
tion, another proposal for enacting a phase gate in this type 
of systems has been put forward 21 . This proposal, while also 
relying on tunneling of Josephson vortices, uses a different 
idea: namely coupling a topological qubit to a flux qubit, and 
effecting a phase gate on the flux qubit instead. 
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